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Abstract. In this paper, we first prove a gradient estimate for the positive 
solutions of the ecumenic degenerate parabolic equation: 

Ut = AF(u), 

writh F'iu) > 0, on a complete Riemannian manifold with Ricci curvature 
lower bound Ric{M) > —k with fc > 0. The second part, we apply the 
gradient estimates to the Fast Diffusion Equations (FDE) and Porous Media 
Equations (PME): 

Ut = A{uP), p > 0, 
to obtain the gradient estimates in a larger range of p than the range of p 
for Harnack inequalities and Cauchy problems in the literature, and also prove 
some Liouville-type theorems for positive global solutions on noncompact com- 
plete manifolds with nonnegative Ricci curvature for the Fast Diffusion Equa- 
tions (FDE) and the Porous Media Equations (PME). 



1. Introduction 

In this paper, we consider a localized Hamilton-type gradient estimates for the 
positive solutions of the ecumenic degenerate parabolic equation 

(1.1) Ut = AFiu) 

on a complete Riemannian manifold (M", g) of dimension n > 1 with Ric{M) > —k 
for some fc > 0. Here F G C^(0, oo), F' > 0, and A is the Laplace-Beltrami 
operator of the metric g. There is a lot of literature on this kind of topics. For 
example, related problems such as the Fast Diffusion Equations and the Porous 
Media Equations have been considered by D.G. Aronson [1], G. Auchmuty and D. 
Bao [2], M.A. Herrero and M. Pierre [6] and S.T. Yau [10]. 

It is well known that, in the study of geometric analysis as well as other elliptic or 
parabolic equations, the gradient estimate and the Harnack inequality play a most 
important role. To begin with, let us review some main results on the gradient 
estimate and the Harnack inequality. The first one is the Harnack-type differential 
inequality for the heat equation by P. Li and S.T. Yau [7], 

Theorem A (P. Li and S.T. Yau [7]). Let Af" be a complete manifold of 
dimension n > 2 with Ricci{M") > — fc for some k > 0. Suppose that u is any 
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positive solution to the heat equation in B(xq,R) x [io — T,to]. Then for a> 1, 



in B{xo,R/2) x [to — T/2,to]- Here Cn depends only on the dimension n and a. 

For the heat equation on compact manifolds without boundary, there is the 
another type gradient estimate by R. Hamilton [5]: 

Theorem B (R. Hamilton [5]). Let M" be a compact manifold without bound- 
ary and with Ricci{M) > —k for some k > 0. Let u be a sm,ooth positive solution 
of the heat equation with u< M for all {x,t) G M" x (0,oo). Then 



The Hamilton-type gradient estimate takes up a significant position in the study 
of the heat equation. However, the classical Hamiltons estimate is a global result 
which requires the heat equation should be posed on compact manifold without 
boundary. Recently, a localized Hamilton type gradient estimate was proved by 
P. Souplet and Q.S. Zhang [9], which can be viewed as a combination of Li-Yau's 
differential Harnack inequality [7] and Hamilton's gradient estimate [5], 

Theorem C (P. Souplet and Q.S. Zhang [9]). Let M" be a complete Rie- 
mannian manifold with dimension n > 1, Ric{M) > —k, k > 0. Suppose that 
F €E C^(0,oo) with F' > 0, and u is any positive solution of the degenerate para- 
bolic equation (1.1) in Qr^t = B{xo,R) x [to ~ T,to] C M x (—00,00). Suppose 
also that u< M in Qr,t- Then there exists a dimensional constant C such that 



for all {x,t) in Qr/2,t/2- 

Moreover, if M" has nonnegative Ricci curvature and u is any positive solution 
of the heat equation on M x (0,oo), then there exist dimensional constants Ci,C2 
such that 



for all x^W andt> 0. 

The first Harnack-type differential inequality dealing with the degenerate para- 
bolic equation (1.1) was attributed to S.T. Yau [10], while the Harnack-type differ- 
ential inequalities have been widely used for Porous Media Equations, see [1], [2], 
[6]. We state Yau's main result as below. 

Theorem D (S.T. Yau [10]). Let M" be a compact Riemannian manifold 
without boundary, Ricci{M) > 0. Suppose that F e C^(0,oo) with F' > 0, c{t) e 
C^(0,oo), and u is any positive solution of the degenerate parabolic equation 



on M". Let a ^ Q be an arbitrary constant. Define a function G on (0, 00) by 

G'{s) =F'{s)/s, and we abbreviate G = G{u), F'' = F''{u), «; = 0,1,2. 
// the conditions below are satisfied: 
(A). |VG|2 - aGt - c{t) <0 att = 0; 







ut = AF{u) 
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(B). (nonlinear condition) the following quadratic inequality holds true for all 
x>0 

l-af ., 2(1 - a) .\ 2 /4(l-a) uF"\ 
> — :^\auF"-^ Lf'\x^+\^ — ^-——\c{t)x 

( 2 uF"\ c^{t) 

then we have for allt>0 that 

|VGp - aGt - c{t) < 0. 

It is natural to seek a localized Hamilton-type gradient estimate for the general 
parabolic equation (1.1) as P. Souplct and Q.S. Zhang [9] did for the heat equation 
on complete manifolds. Under some strong assumptions a localized Hamilton-type 
gradient estimate for the equation (1.1) was obtained by L. Ma, L. Zhao and X. 
Song [8]. The first main result of this paper is to generalize the rc;sult in [8] to 
the following Hamilton- type gradient estimate for the equation (1.1) under some 
weaker assumptions: 

Theorem 1.1. (Gradient Estimate). Let M" be a complete Riemannian man- 
ifold with dimension n> \, Ric(M) > —k, k > 0. Suppose that F G C^(0, oo) with 
F' > 0, and u is any positive solution of the degenerate parabolic equation (1-1) in 
Qr^t = B{xo,R) X [to — T, to] C M X (—00,00). Define a function G on (0, 00) 
by G'{s) — F'{s)/s. We denote by U <Z (0,oo) the value range of u. If there exist 
nonnegative constants K, a, (5, r and 7 satisfying: 

(A) . F'{s)<K, WseU; 

(B) . a - G{s) >d>0,y s&U; 

(C) . (Nonlinear Condition) 

2+^(2-(--l)^^)>7>0, yseU 



then there exists a constant C{n, S, K,T,"f) depending only on n, S, K, r and 7 
such that\/{x,t) G Qr/2.t/2> there holds 

When n = 1, the Ricci curvature lower bound k vanishes. 

Remark 1.1. One might compare the above theorem with Theorem 7 in [8]. the 
nonlinear condition ( C) here is much weaker than the one in [8] . Our Nonlinear 

Condition (C) can apply to a much larger class of equations. One may see it from 
those results on the Fast Diffusion Equations and the Porous Media Equations in 
the second part of the paper. 
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The second part of this paper, we study the Fast Diffusion Equations and tlic 
Porous Media Equations: 



on a complete Ricmannian manifold (M" , g) . Firstly for p < 1 , wc generalize the 
gradient estimates for (1.3) in [8], where the theorem holds only for dimension 
n = 2 or 3 with much shorter range of p (only for p € {1 — l/^/n, 1) for n = 2, 3 in 
[8]), as the following result: 

Theorem 1.2. Let M" be a complete Riemannian manifold with dimension n>l, 
Ric{M) > —k, k > 0. Suppose that u < M is a positive solution of the Porous 
Media Equation (1.3) in Qr^t = B{xo, R) x [to — T,to] C M x (— oo, oo), where 



for all {x,t) in Qr/2,t/2- When n = 1, the Ricci curvature lower bound k 

vanishes. 

Remark 1.2. One should notice that the range of p here is (1 — ;;j^,l), while 
previous results for the Fast Diffusion Equations on R" requires p€((l — ^) + ,l), 
for example, a survey by D.G. Aronson [1], the Cauchy problem, by M.A. Herrero 
and M. Pierre [6], the Harnack-type inequalities by G. Auchmuty and D. Bao [2]. 
We can see that for n > 3, the range ofp for our gradient estimate is larger than the 
range ofp in previous results [1], [2], [6]. [8]. Our gradient estimate will be a useful 
tool to study related problem for the Fast Diffusion Equations in this large range 
of p, with which one couldn't deal before even on R". We shall study the Cauchy 
problem, the Harnack-type inequalities for the Fast Diffusion Equations (p <1) in 
another paper. 

On a noncompact manifold with nonncgative Ricci curvature, an immediate 
application of Theorem 1.2 is the following time-dependent Liouville theorem of the 
Fast Diffusion Equations, which generalizes Yau's celebrated Liouville Theorem for 

positive harmonic functions. 

Theorem 1.3. (Liouville theorem) Let be a complete, noncompact manifold 
with nonnegative Ricci curvature. Let u be a positive ancient solution, a solution 
defined in all space and negative time, of the Porous Medium Equation for 1 — < 
p < 1 , and L{s) £ C(R) be any strick increasing function with L{s) oo as 
s ^ oo, such that 



near infinity. Then u is a constant. 

Remark 1.3. One might see that the growth condition in the spatial direction in 
Theorem is very weak, since we might choose 

L{s) = exp(exp(- • • (exp(s)) • • • )) 



(1.3) 



1 < p < 1, for n > 1. 

n + 6 

Then there exists a constant C (n, p) depending only on n and p such that 



(1.4) 
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with I exp for any I > 0. Note that one might write any positive harmonic function 
v{x) as a positive global solution of lS.{vP) = with u{x) — v{xy/P, hence Yau's 
celebrated Liouville Theorem for positive harmonic functions is a special case of 
Theorem 1.3 for time dependent positive solutions of the Fast Diffusion Equations, 
while one couldn't do this for the Heat Equations(see examples in [9]). 

Secondly, we shall study the Porous Medium Equation for p > 1. We have the 
following gradient estimate for n = 1: 

Theorem 1.4. Suppose that u is a positive solution of the Porous Media Equation 
(1.3) in Qr.t = B{xo,R) x [to — T,to] C M x (— cxd,oo), with range(w) = [m,M]. 
Let G{u) = i^vP-^, a = :^Mp-\1 + 6) with some constant 5 > Q. Then for 
any p> 1, there exists a constant C{p) depending only on p such that 

|V,G(«(.r,i))| <c{p)(^^^ ' ^ ^ 



a-G{u{x,t)) ^'\y5R VMp-'^ST J ' 
for all {x,t) in Qr/2,t/2- 

An immediate application of this theorem is the following time-dependent Li- 
ouville theorem for the Porous Medium Equation with p > 1 and the Liouville 
theorem for positive solutions of A(uP) = when n — 1^ 

Theorem 1.5. (Liouville theorem) Let u be a positive ancient solution, a solu- 
tion defined in all space and negative time, to the Porous Medium Equation (p> 1) 
on R, such that 

u{x,t) = o(d{xy/^p-^^ + 

near infinity. Then u is a constant. 

And for dimM > 2, we have the following gradient estimate: 

Theorem 1.6. Let M be a complete Riemannian manifold with dimension n >2, 
Ric{M) > —k, k > 0. Suppose that u is a positive solution of the Porous Media 
Equation (1.3) in Qn,T = B{xo,R) x [to — T,to] C Mx (—00,00), wii/i range(u) = 
[m,M]. Let G{u) = ^^u'^~^, a. = ■^^MP~^{1 + 5) with some small constant 
< (5 < // the following pinch condition on m, M holds 

Ap \ 




l + 5\{n-l){p-l) 7' 
then there exists a constant G(n,p) depending only on n and p, a 

(n-l)(p-l) MP-^n + d) -mP-^ „ 

7 = 2p-^ '-^ ^—3^ >0, 

2 mP ^ 



such that 

\Va,G{u{x,t))\ 



<G{n,p)\ 



a-G{u{x,t)) ' 'y^ySR ^^5Mp-'^T 

for all {x,t) in Qr/2,t/2- 
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The rest of this paper is organized as foUows. In section 2, we prove Theorem 
1.1 following the argument in [8], which can be regarded as a combination of the 
analysis in [7] and [9]. In section 3, we apply the Theorem 1.1 to study the Fast 
Diffusion Equation and the Porous Media Equations and prove Theorem 1.2-1.6. 
Here and later B.icci{M) is the Ricci curvature and a manifold is complete if every 
geodesic extends to infinity. 

We shall use the Einstein summation convention for indices i, j, k, etc. In 
particular, we shall use the short-hand notation gfj = gijgij = X^"j=i 9ij- 

2. Hamilton-type gradient estimates 

In this section, we shall prove Theorem 1.1 and show that the Harnack quantity 
trick introduced in the fundamental work of P. Li and S.T. Yau [7] can be adapted 
to study the Hamilton-type estimate. Our argument can be seen as an improvement 
of the argument in [8], which can be regarded as a combination of the analysis in 
[7] and [9] . We first show a technical linear algebra lemma on symmetric matrix. 

Lemma 2.1. Let A = (uij) be a nonzero nxn symmetric matrix with eigenvalues 
{Aft}, for any a, 6 € R, one has the following properties: 

n n 

(a). \A\'=Y,al=tv{AA'') = J2>^h 

i,j=l fc=l 

n 

(6). max (aA + htr AIn) {v, v) = aXi + b A/^, for some 1 <i <n 

n 

min (oj4 -|- 6trA7„) {v, v) = aXj + b\2^ky for some 1 < j <n. 

f\ \aA + btrAIn, ^l^ , i,n2 , / .s,2 

• . e^?^ , , TTi ^' ^ ={a + b) +{n- 16 . 

AeS{n);\v\ = l I \A\ J 

Proof, (a) follows from direct computation and A symmetry. 

(b) follows from the facts that max|i,|^i (^aA + btrAIn) {v, v) and min|„|^i (oA + 
btTAIn){v,v) arc the maximal and minimal eigenvalue of aA + 6trA/„, and the 
eigenvalues of aA + btiAIn are {aXi + bJ2k=i •^ft}?=i- 

To prove (c), apply (a) and (b), we have 



max 

AeS{n);\v\ = l 



aA + b\vAI„, 

\v,v) 



max [aAi + Afe]^ 



By Lagrangian multiplier method in Calculus, the extremums of /(Ai, • • • , A„) = 
aAi + Afe under constrain J2k=i ^k~^ 

n 

- [{a + bf + (n - l)62]i/2 < aAi + 6^ Afe < [(a + bf + (n - 1)62]V2. 

fe=i 

□ 

The rest of this section, we shall prove Theorem 1.1. We shall define a quantity 

w{x, t) and the starting part of the argument will be to deduce a differential in- 
equality on w{x, t), then use the well-known cut-off function of Li- Yau [7], to derive 
the desired bounds. 



GRADIENT ESTIMATES AND LIOUVILLE THEOREMS 



7 



Let (j) — liiu. Since u is a solution to the equation ut — A{F{u)), simple 
calculation shows 

(2.1) = A(G(w)) + VG{u) ■ Wcj). 

Writing g{(])) = G{e'^), and multiplying (2.1) by g'{(p)a,nd some elementary compu- 
tations, we then get 

gt = g'Ag+\Wg\^, Vg = g'y<t>. 
g'{cl>) = G'(e*)e* = F'(u). 
g"[(f>) = F"{e^)e't' = F"{u)u. 

Set 

(2.2) u; = ^a;,t) = |vin(a-5)|' = ^^^^, 

(a - 9) 

and we first derive a differential inequality for w, to which we apply the maximum 
principle. 

Lemma 2.2. w satisfies the following differential inequality: 

g'Aw — wt > Lw^ — 2g'kw — Li\/g ■ Vw, 

where L and Li are some functions give in (2.10) and (2.11). 

Proof. After some elementary computations in local orthonormal system as in 
[8], we get that 

= + 2-1^^1'^* 



(a - gf {a - g)^ 
^ Vg-V(.9-A.g+|VflH ^ JVg\^g'Ag+\Vg\^) 

{a - gf {a - gf 

2g'Vg ■ VAg 2AgVg ■ V g' 2Vg ■ V|Vgp 2g'|VgpAg 2|V5|^ 



{oi-gY {(x-gf {oi-gY {a-gf {a-gf 

(2.3)= 2g' , + 2^^ '^■'^''^^ + 4 ^^^'-'^^ + 2g' '^^''^f + 2 '^^1' 



(a - gY g' (a - gY - gY {a - gY (q - gY 

(2.4) = 



{a-g^Jj {a-gY {a-gY 



Aw = Wii — 2 



\ , r,( gfgj 



" \{a-gY)i \{a-gf 
(2.5) = 2 I 2 I g ^'^'^"^^ I 2 *^^*^^^ I 6 



(a-ff)^ (a-ff)^ (a -5)^ (" - ff)^ {ot-gY 
By (2.3) and (2.5), we obtain that 

g'Aw-Wt = 2g'^^^ + 2g'^l^l^-^l§^ + 8g' 



{a - gY (a - gY [a - gf 

gigm , . , |V.9|^ J' \^g?Ag ^ \^g\^ 



^ (a - 5)2 + (a - g)4 2-^, _ 2-^^ _ 

Bochners identity implies that 

gigin - gjgiij = gjigju - gm) = Rugigj = Ric{Vg, Vg) 
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where Rij is the Ricci curvature tensor. Therefore we arrive at 



(2.6) 



= 2g' 



+&g 



{oc-gf 



+ 2g' 



,Ric{Wg,Vg) , gigijgj 9i9ij9j 



(a-gy 



,g" \Vg\'Ag ^ \Vg\^ 



{o^-gT 9' {oi-9Y {oi-9Y 



[a -9? 



Recalling (2.4), we have 



(2.7) 



(a - gf {a - gf 



Adding ^2 — — r}^^ x(2.7) with (2.6), where is a parameter function to be 
determined later, we conclude that 



g'Aw — Wt 



V- 



9t 



{a - ff)2 
+ (25' 



+ 2 



2.g' 



a — g 



g"\grgm ^g" iVg^Ag 



2-2,^ 
9 



g' J {a- g)2 g' {a - g)^ 
|V5|^ 



(2.8) 



+2g' 



, Ric{Vg,Vg) 
{a - gY 



{a-g))^ 

J {a- gY 

V a-g 'a'J ^ 



Denote / = 2g'/{a — g), b = g"/g', A = [gij) and v = V^/IVflf], then we have 
(/fj = \A\^, 9i9ij9j = A{v,v)\Vg\^ and Ag = tiA. Prom definition of w, the right 
side of (2.8) can be written as 



> 



(2.9) 



Vr^^ + 2(a-5) 

{a - gr 



{f-vb) 



A{v,v) tvA 



1^1 

a-g 



-w 



+(" - 5')(/ + 2 - 2r]bjw'^ + 2g'Ric{v,v)w - {2 - f - r]b^Vg ■ Vw 



25' 



1^1 



{oi-g) 



A(v,v) ArA\ 



H 2 



+(a - <?) + 2 - 2,76 - i 



{f-rih) 



A{v,v) trA 



w 



+2g'Ric{v, v)w- (2- f- 2r]b^\7g ■ Vu; 



2-2r]b 



(/ - Vb 



A{v,v) tiA 



w 



+2g'Ric{v, v)w- (2- f - r]b^ Vc/ • Vw 
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Apply Lemma 2.1 to (2.9), we have 



g'Aw-Wt > ^(^f+(2-2rjb)f-[f-T,b-by-{n-l)b^^ 

+2g'Ric{v, v)w- ^2 - / - 'qb^ Vg ■ Vw 

= ^{2il + b)f-{n-l + [rj + ir)b')w' 

+2g'Ric{v, v)w- ^2 - / - rjbj Vg ■ Vw 

Next we estimate the coefBcient L oi xiP . Pick 77 = —1, We can bound L as the 
following: 

(2.10) L={a-g) (2(1 + 6) - (n - l)y) > (a - 5)7 > 0, 

since 2(1 + &) — (n — l)y > 7 > from our Nonlinear Condition (C) in Theorem 
1.1. We estimate the coefficient Li of Vg • Ww as 



2-f + b 



<2 + r + f. 



(2.11) = 
Hence from Ric(M) > —k, we have 

(2.12) g'Aw -wt> Lv? - 2g'kw - LiVg ■ Vw 

□ 

Now we can apply maximum principle on the differential inequality (2.12) to 
prove our gradient estimate (1.2). We will follow [8] and [9] to use the well-known 
cut-off function by Li and Yau [7] to show Theorem 1.1. We caution the reader that 
the calculation is not the same as that in [7] due to the difference of the first-order 
term. 

Proof of Theorem 1.1. Let "if = "^{x, t) be a smooth cut-off function supported 
in Qr,t, satisfying the following properties: 

(1) . * = ^{d{x,xo),t), * = 1 in Qr/2,t/2, < ^' < 1; 

(2) . ^ is decreasing as a radial function in the spatial variables; 

(3) . %^<%, |^<§,whenO<a<l; 

(4) . ^ < ^ 

Then, from (2.12) and a straightforward calculation, one has 

V\E' 

(2.13) g'A{^w) - (*w)t + L1V5 • V(*w) - 2g'— ■ V(*«;) 

> L^w^ - 2kg'^w + g'A^w - ^tw + Li{Vg ■ W^)w - 2g' ^ ^' w. 

Hereafter we always use C to denote different constant depending on n only. We 
obtain the upper bounds for each term of the right-hand side of (2.13) as did in P. 
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Souplet and Q.S. Zhang [9] 



(2.14) 
(2.15) 
(2.16) 



2kg'^w 



2g'^-—^w 



< 



L*w;2 H 1 < -L*u;2 + CA;^^ 



/2 



6 



2 , W 

2L* 

„'2 



1 



C 1 



1 r T 2 65 |V*|4 1 . . 2 5 



Li(V5- V*)m; 



Li(V5f V*)-:^:r-^W 



^^,,,3/2 

IV.9I 



(2.17) 



^ Irr 2 CLt(a - ff)4 |V*|4 

6 *3 
^ 1,, 2 L\{a-gY 



Here we make use of Young's inequality, 



ab< \ , 

p q 



Vp, (7 > with — I — = 1. 



Furthermore, by the properties of ^ and the assumption on the Ricci curvature, 
one has 

l-^'A^-w;] < 5'[92^ + (n-l)^ + 9^*a^ln(Vg) 



< 5'[|a2*|+2(n-l)M 

< -L*?«^ + — 
-6 L 



+ 4(n- 1)2^^1^ + fcl 



(2.18) 



< ^L'i/'w^ + C 
6 



1 1 



Inserting (2.14)-(2.18) into the right-hand side of (2.13), we deduce that 
g'A{^w) - {^w)t + Z/iVfif • V(*u;) - ^9 ■ V(*w) 



- 6 



(2.19) 



L 

> — 
- 6 



Li{a-gf 1 g'^ 1 11 
L\{a-gY 1 



/2 



.'2 1 



L4 



^ ^ 5^ 2 



RecaUing that i > 27(a - 5) > 0, Li < 2 + r + / and / = g' /{a - g), that is 



g'A{-^w) - {■^w)t + LiVg ■ V(*m;) - 2.g'— • V(*w) 



> — 

- 6 



(2.20) > 



-Pw^ - C(7) I 



(l + / + r)^ 
' (l + / + r)^ 

ii4 



+ 



+ 



9 1 /2,2 



(a - 5)2 T2 



2;„2 
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Suppose that the maximum of (^w) is reached at {xi,ti) G Qr.t- By [7], we 
can assume, without loss of generahty, that xi is not in the cut-locus of M . Then 
at this point, one has A(*ui) < 0, (*u')t > and V(*w) = 0. Recalling that 
ct — g>5,Q<g'<K and / < K/5, we have 

By assumption, the maximum of {^w) is reached at {xi,t\) G Qr,t, which 
implies that for any (a;, t) G Qr,t 

(*«.y(a:,i) < {^wy{xi,ti)< (^w^yxuh) 

Noticing that ^{x,t) = 1 in Qji/2,T/2 and w = |Vfif|^/(a — g)^, we finally have 
proved 



(a-g)2 - S^R^ ST S 

which is exact the conclusion of Theorem 1.1. □ 



3. Some applications 

In this section, we shall study the Heat Equations, the Fast Diffusion Equa- 
tions and the Porous Media Equations, and obtain some Hamilton-type gradient 
estimates by applying our Theorem 1.1, then prove some time-dependent Liouville 
Theorems on noncompact complete manifolds with nonnegative Ricci curvature. 
As a corollary, we obtain Yau's celebrated Liouville theorem for positive harmonic 
functions: any positive harmonic function on a noncompact manifold with nonneg- 
ative Ricci curvature is a constant function. 



3.1. Heat Equations: Let M" be a complete Riemannian manifold with dimen- 
sion n > 1, Ric{M) > —k, k>Q. Suppose that w < M is a positive solution of the 
heat equation 

ui = Au 

in Qn.T = B{xq,R) x [to — T.ta] C M x (— oo.oo). We may look as the case of 
F(u) = u in Theorem 1.1. Choose G{s) = Ins and let a = 1 -|- InM, K — 1, t — 1 
and 7 = 2 in Theorem 1.1, we obtain the Theorem C by P. Souplet and Q.S. 
Zhang in [9]) as a direct corollary of our Theorem 1.1. 

Corollary 3.1. (Theorem 1.1 in [9]) Let M" be a com,plete Riemannian manifold 
with dimension n > 1, Ric(M) > —k, k > 0. Suppose that F e C^(0,cx)) with 
F' > 0, and u is any positive solution of the degenerate parabolic equation (1-1) in 
Qb.,t = B{xo, R) X [to — T, to] C M x (— oo, oo). Suppose also that u< M in Qr,t- 
Then there exists a dimensional constant C such that 



u{x,t) ~ \R \/T 
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for all {x,t) in Qr/2,t/2- 

Moreover, if M" has nonnegative Ricci curvature and u is any positive solution 
of the heat equation on M x (0,C!o), then there exist dimensional constants C\,C2 
such that 

\vM.M^^i(^^uix^\_ 

U[X,t) ^/ty U{x,t) J 

for all X e M" and t > 0. 

Using the above gradient estimates for positive solutions of the Heat Equation, 
P. Souplet and Q.S. Zhang in [9]) proved the following time-dependent Liouville 
Theorem: 

Corollary 3.2. (Theorem 1.2 in [9]) Let M be a complete, noncompact manifold 
with nonnegative Ricci curvature. Then the following conclusions hold. 

(a) Let u be a positive ancient solution to the heat equation (that is, a solution 

defined in all space and negative time) such thatu{x,t) = e°(''(^^+\/i*l) near infinity. 
Then u is a constant. 

(b) Let u be an ancient solution to the heat equation such that u{x, t) = o{[d{x) + 
near infinity. Then u is a constant. 

As discussed in [9], both growth conditions of the above theorem in the spa- 
tial direction are sharp, by the some simple examples. Hence one couldn't obtain 
Yau's celebrated Liouville theorem for positive harmonic functions directly from 
the above time-dependent Liouville Theorem for heat equation. 



3.2. Fast Diffusion Equations: Let M" be a complete Riemannian manifold 
with dimension n > 1, Ric{M) > —k, k > 0. Here we consider the Fast Diffusion 
Equation 

ut = A(iiP), p < 1 

on M X (—00, oo). We have the following localized Hamilton-type gradient estimates 
for the positive solution on M x (—00,00): 

Theorem 3.1. Suppose that u < M is a positive solution of the Fast Diffusion 
Equation 

Ut = A{uP) 

in Qr,t = B{xo, R) x [to - T,to] C M x (-00, 00), where 

(3.1) l-^<p<l, forn>l. 

n + 3 

Then there exists a constant C depending only on n and p such that 

u{x,t) ~ \R ' 

for all {x,t) in Qr/2,t/2- When n = 1, the Ricci curvature lower bound k 
vanishes. 
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Proof. For < p < 1, we have F{s) = s^, and choose G{s) = ^r^-s^ ^, and 
let a = in Theorem 1.1. We have F'{s) = psf-\ F"{s) = p{p - l)sf-2, then 
Condition (C) in Theorem 1.1 becomes 

T >1— p 

Ap- (n- 1)(1 -p) > 0. 
which is equivalent to, 



1 - 



n + 3 



< p < 1. 



Let 7 = ((n + 3)p — (n — 1)) /2 for given p in admission range (3.1). We follow the 
proof of Theorem 1.1, wc have 

9^^yf-\ g'=puP-\ g" =pip-l)uP-^ 
/=^=2(l-p), b=^=p-l, 

L = l^uP-\ L,=3-p. 

We follow the proof of Theorem 1.1 until (2.19) with the modification that the 
constant C (n, p) here depends only on n and p, we insert the above quantities into 
(2.19) to obtain 

2g'A{^w) - {^w)t + LiVg ■ V(*u;) - • V(*«;) 



> 



1 



^w'-C{n,p)\ — + 



1 



By the same maximum argument in proof of Theorem 1.1, we have 

■ 1 M2(1-p) 



T2 



which implies 



<Cin,p)[^ + 



+ k 



Then the conclusion of Theorem 3.1 follows easily from the fact that 

iVffI _ \VG{u)\ _ \Vu\ 
-g -Giu) ^ u 

□ 

An immediate application of the above gradient estimates is the following time- 
dependent Liouville theorem for the Fast Diffusion Equations on a noncompact 
manifold with nonnegative Ricci curvature: 

Theorem 3.2. (Liouville theorem) Lei M" he a complete, noncompact manifold 
with nonnegative Ricci curvature. Let u he a positive ancient solution, a solution 
defined in all space and negative time, of the Fast Diffusion Equation for 1 — < 
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p < 1 , and L{s) G C(R) be any strick increasing function with L{s) — J- oo as 
s — >■ oo, such that 

u{x,t) = o(^L{d{x)) + \t\^/^^-P^^ 
near infinity. Then u is a constant. 

Proof. Since L(s) is a strick increasing function with L{s) — )• oo as s ^ oo, 

there is inverse function H{s) of L(,s) which is also a strick increasing function with 
H{s) — )■ oo as s — >■ 00. Fixing (xo,to) in space-time and using Theorem 3.1 for u 
on the cube Q{^H{Ry(^-P^), R'^) = B{xo, ^H{Ry(^-P^)) x [to - i^^^o], and the 
M in Theorem 3.1 is the maximum value on the double cube 0(iJ(i?2/(i-p)), 2i?2), 
by our assumption on the growth condition of the function u at infinity, 

M^(«va-.)).2H= =o(L(i?(i?2/(i-rf)) + ii!2/(i-rf) =o(i?^/(i-^)). 

Hence by Theorem 3.1, we have 

|VM(a;o,io)| ^ ^/ 1 



(imi?2/(i-P))+:^''^^0' 



u{xo,to) - \\H{RV(.^-p)) 

Letting i? — >■ oo, it follows that |Vu(a;o,to)| = 0. Since {xQ,to) is arbitrary, one sees 
that u = c. □ 
As a corollary of the above time dependent Liouville theorem for Fast Diffu- 
sion Equation, we obtain Yau's celebrated Liouville theorem for positive harmonic 
functions on a complete, noncompact manifold with nonnegative Ricci curvature: 

Corollary 3.3. (Yau's Liouville theorem for positive harmonic functions) 

any positive harmonic function on a noncompact manifold with nonnegative Ricci 
curvature is a constant function. 

Proof. The proof follows immediately from the above time dependent Liouville 
theorem. Let t; be a positive harmonic function. Choose an p with 1 — <p <1, 

then u{x) = v{xY^p is a positive solution of A(u^') = 0, which can be regarded as 
a time independent positive solution of the corresponding Fast Diffusion Equation. 
Define 

L{s) = 8 max u{x) + s 

d{x)<s 

It easy to see L{s) is a strick increasing function with L{s) — > oo as s — >■ oo, and 
u{x, t) = o{L{d{x))) near infinity. Follow from the above Theorem, we have u must 
be a constant, which implies v must be a constant. □ 

3.3. Porous Media Equations: Let M" be a complete Riemannian manifold 
with dimension n > 1, Ric{M) > —k, k > 0. Here we consider the Porous Media 
Equations 

ut = A{uP), p > 1 

on M X (— oo, oo). 

Firstly for dimension n = 1, we have the following localized Hamilton- type gra- 
dient estimates for the positive solution on M x (— oo, oo): 

Theorem 3.3. Suppose that u is a positive solution of the Porous Media Equation 
on M" dimension n = l 

Ut = A(nf) 



GRADIENT ESTIMATES AND LIOUVILLE THEOREMS 



15 



in Qr,t = B{xo,R) x [to — T,to] C M x (—00,00), with range('«) = [m,M]. Let 
G{u) = ct = ^^M^'~^(l + 15) with some constant ^ > 0. Then for any 

p> 1, there exists a constant C{p) depending only on p such that 

\W,G{u{x,t))\ + 5 ^ 1 ^ 



a-G{u{x,t)) ^'ySR y/MP-'^6T J ' 

for all {x,t) in Qr/2,t/2- 

Proof. Wc have F{s) = and choose G(.s) = j^-''^^^ in Theorem 1.1. We 
have F'{s) — ps^^^, F"{s) = p{p — l)s^'^^,and K — pMP^^, and a, as defined in 
above, then Condition (A) and (B) are satisfied and Condition (C) in Theorem 1.1 
becomes 

T>p-1, 

2p > 7 > 0, V s e [m, M] 

Condition (C) is satisfied if we choose 7 = 2p and t = p—1. As proof of Theorem 
3.1, we have 

9 = g' = puP-\ g" = p{p - 1)uP-^ 

/ ~ a-g — "^yy ■LJmp-i(1+5)-uP-i - 5 

, — aZ 



b=^=p-l 



L = 'y{a-g) = ^ Mp-\1 + 6)- u^-' 



> 2p_^p-lj 

— i-p 



L,=i+p+f<i+p+2{p-i) M.-i(ir/)-«p-^ <i+p+ 

By the same argument in proof of Theorem 3.1, we have 

< c{p) f i + T^ol < cip) f ii±i^ + ' ^ 



which implies 



(a - gf - \ B? MP-'^ST J 
Then the conclusion of Theorem 3.2 follows easily from the fact that 

|Vg| _ IVGWI 

a — g a — G{u) 

□ 

An immediate apphcation of this theorem is the following time-dependent Liou- 
ville theorem for the Porous Medium Equations on R, 

Theorem 3.4. (Liouville theorem) Let u be a positive ancient solution, a solu- 
tion defined in all space and negative time, to the Porous Medium Equation (p> 1) 
on R, such that 



i{x,t) = o(d{xy/^p-^^ + 
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near infinity. Then u is a const,ant. 

Proof. By our assumption the function u satisfies u{x,t) = o{d{xy^''P~^^ + 
|i|i/(p-i)) near infinity. Fixing {xo,to) in space-time and using Theorem 3.3 for u 
on the cube Q{R, R) = B{xo, R) x [to — R, to], we have 



< C{p) ( -TTT + ) M"-^ < C{p, S) ( 



u(xo,ta) - '-^'V SR VMp-^ST^ " ' V i? V T 

where M is the maximum value on the double cube Q(2R, 2R), by our assumption 
on the growth condition of the function u at infinity, 

Hence we have 



Letting i? — )• oo, it follows that |Vu(a;o,to)| = 0. Since {xo,to) is arbitrary, one sees 
that u = c. □ 
Secondly for dimension n > 2, we have the following localized Hamilton- type 
gradient estimates for the positive solution on M x (—00,00): 

Theorem 3.5. Let M be a complete Riemannian manifold with dimension n>2, 
Ric{M) > —k, k > 0. Suppose that u is a positive solution of the Porous Media 
Equation 

ut = A{uP) 

in Qr,t = B{xo,R) X [to — T,to] C M x (—00,00), with range(w) = [m, M]. Let 
G{u) = ^uP-\ a = ^MP-i(l + S) with some small constant 0<S< // 
the following pinch condition on m, M holds 

then there exists a constant C(n.p) depending only on n and p, and 



^ {n-l)(p-l)MP-\l + 5)-mP-^ ^ 



mP- 



such that 




a — G{u{x, t)) 
for all {x,t) in Qr/2,t/2- 

Proof. We have F{s) = sP and choose G{s) = ^^sP~^ in Theorem 1.1. We 

have F'{s) = psP~^, F"{s) = p{p — l)sP~^,and K = pMP~^, and a, as defined in 
above, then Condition (A) and (B) are satisfied and Condition (C) in Theorem 1.1 
becomes 

T>p-1, 

2p - (n-DiP-i) M^-^(i+6)-.^-^ ^ ^ ^ 0^ Vse[m,M] 
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which is equivalent to, 

T > p- 1, 

Let T = p — 1 and 7 = 2p — (iiilLKp^:!! M — ^^^i"™ — > 0, the above condition as 
< (5 < is equivalent to our pinch condition (3.2). 
As proof of Theorem 3.3, we have 

9 = g' = puP-\ g" = pip - 

J — a-g — "^yy -LJmp-i(1+5)-uP-i - 5 



b=^=p-l 



L = j{a-g) = ^ MP-i(l + 6)- u^-^ > 



Li = 1 + p + / < 1 + p + 2(p - 1) < 1 + 

By the same argument in proof of Theorem 3.3, we have 

( [5 + lY 1 A;^' 



2(P-1) 



which implies 



(a -5)2- ' ''^'K^-^d'^E? -ySMP-'^T 
Then the conclusion of Theorem 3.5 follows easily from the fact that 

|Vg| _ |vg(»)l 

a — 5 a — G{u) 

□ 
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